Quantum motion of charges in 7r-conjugated polymers: kinkiness induces localization 
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We develop a model for charge trapping on conjugated polymer chains using a continuous repre- 
sentation of the polymer. By constraining the motion to along the chain, we find that kinks along 
the chain serve as points of attraction and can induce localization and spontaneous buckling of the 
chain. We implement this in a model system of a conjugated polymer donor and fullerene acceptor 
to mimic a polymer based photovoltaic system. We find that geometric fluctuations of the poly- 
mer give rise to deep potential energy traps at kinks which are often stronger than the Coulomb 
interaction between the electron and hole. 



One of the persistent questions arising in the modeling 
of polymer-based solar cells is how the morphology and 
dynamics of the polymeric molecules comprising the cell 
relates to overall efficiency of the cell to separate photo- 
excitations into mobile charges. On the short time-scale, 
we can consider the charge-separation event in a frame- 
work of a fixed nuclear geometry. For certain, the effi- 
ciency and rate of exciton fission is intimately dependent 
upon how donor and acceptor species are in contact with 
each other. Once exciton fission has occurred to produce 
a bound charge-transfer pair, one needs to consider how 
the fluctuations and dynamics of polymer chains them- 
selves play into further separating this state into mobile 
charge carriers that can be collected. However, in or- 
ganic conjugated polymer semiconductors, the idea of a 
"mobile charge" is somewhat of a misnomer since disor- 
der, strong-electron phonon coupling, and chemical de- 
fects can lead to charge trapping and drastically diminish 
mobility. [U [2] 

The underlying questions are whether or not the tt- 
conjugation persistence length stretches over the entire 
length of a polymer chain, and will the mobility of a 
free charge on the chain be limited by geometry of the 
chain itself. To address these, we consider a conjugated 
polymer to be a continuous filament in R 3 described by 
the vector j(s) where s is the arc-length on which a mo- 
bile charge with mass m is bound via a confining po- 
tential. While the polymer may be bent and kinked, we 
assume that the tt electronic system is continuous along 
the chain and that the particle can freely move from one 
end to the other. In the limit that the polymer is fully 
extended, the free-charge would behave as a particle in a 
one-dimensional box. 

However, when introducing bends and kinks in the 
polymer, one needs to consider the effect of constrain- 
ing the motion to follow the three-dimensional geometry 
of the polymer itself. In general, the quantum motion 
of the charge satisfies the time-dependent Schrodinger 
equation iHip(q) = Hip(q) with H given by the kinetic 
energy operator in 3-dimensions. Da Costa showed that 
these conditions lead to a reduced Schrodinger equation 
for the longitudinal motion of a quantum particle along 



j(s). [3j [4] We briefly summarize this result here. To 
constrain the motion to along a filament, we define a 
coordinate frame in terms of the tangent to the curve, 
t(s) = 7'(s), and vectors fi2 and hs. These unit vectors 
define a local orthonormal coordinate frame centered at 
7(5) and are related by rotation, T(0) about i to the 
normal and binormal vectors n and b. Any point in the 
ri2, fis plane can be written as R(s) = j(s) + #2^2 + #3^3- 
Vectors t, n, and b satisfy the Frenet-Serret equations 
along j(s): J f (s) = F-'y(s) where F is the anti-symmetric 
Frenet-Serret matrix with depends upon the curvature, 
k and torsion r. 

Note that the derivative of the rotation angle with re- 
spect to the arc- length is equal to the torsion, 0' = r. 
Using R(s), one then finds the covariant metric ten- 
sor is diagonal with elements G\\ = (1 — ft/) 2 , and 
G22 = G33 = 1 where / = ^2 cos 0(s) + q3s'm6(s) and 
0(s) is the rotation of the frame about t. Including 
this rotation of the frame relative to the Frenet-Serret 
frame brings the covariant metric into a diagonal form 
at all points along the curve. Inverting G^ v and writ- 
ing the kinetic energy operator in terms of the Laplace- 
Beltrami operator, allows one to construct an effective 
one-dimensional Schrodinger equation for the free, but 
constrained, motion along the curve 

The original goal of da Costa's work was to understand 
the operator ordering problem when solving quantum 
problems in curved metrics; however, for the purposes 
of understanding charge localization and motion in con- 
jugated polymer-based devices, it provides a useful way 
to connect the influence of the geometry of the the chain 
with the dynamics of a mobile charge along the chain. 

Consider the case of a single quantum particle bound 
to a flexible rod given by, 7(5) with bending modulus, A. 
We can write the total ground-state energy is a functional 
of k(s) and the longitudinal wave function, ip(s) 

mi>] = ^ £ K 2 (s)ds + W\HhM, (2) 
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FIG. 1: Ground-state energy (a) and probability density (b) 
for a quantum particle on a flexible parabolic rod with bend- 
ing modulus A. Our units are such that E Q — h 2 7r 2 /2mL 2 = 1 
is the ground state energy of a particle in a ID box. 



subject to the constraint that the length of the rod is 
fixed L = J Q L \y(s)\ds and ip is normalized to unity with 
boundary conditions ^(0) = i/)(L) = 0. Bending the 
rod about the center increases the curvature and lowers 
energy of the quantum particle bound to the rod but 
increases the strain energy. At some point, the two come 
into balance and we can minimize the total energy self- 
consistently. As a trial curve, we take a parabolic curve 
7(s) = {a(s — 1/2), b(s — 1/2) 2 , 0} with a and b chosen so 
that 7(5) is a fixed length and minimize the total energy 
with respect to ip and the curvature. In Fig. |TJl we show 
energy vs. bending modulus for such a system. The 
dashed line corresponds to the ID box limit where the rod 
is perfectly straight. As the bending modulus decreases, 
the rod is increasingly bent about its middle and the 
quantum state becomes increasingly localized as seen in 
Fig.[lJ>. 

This flexing and localization is spontaneous for all val- 
ues of the bending modulus and should be a universal 
property of all quasi-one-dimensional quantum systems. 
That stated, 7r-conjugated molecular systems and car- 
bon nanotubes are very stiff; however, it may be pos- 
sible to see a spontaneous buckling effect in sinusoidal 
or S'-shaped systems. For purposes of this paper, we 
are focused upon how geometric fluctuations in a con- 
jugated polymer chain can contribute to quantum local- 
ization and charge separation in optical-electronic and 
photovoltaic materials composed of conjugated polymer 
molecules. 

We developed a straightforward model for simulat- 
ing the dynamics in a model bulk hetero junction photo- 
voltaic system consisting of Cqo derivatives blended with 
poly-thiophene chains. Such systems and their deriva- 
tives are of considerable interest for fabricating highly- 
efficient photovoltaic devices based upon organic ma- 
terials. [5-9 Since we are primarily interested in the 
coarse-grained as opposed to atomistic motion of the sys- 
tem, we treated the Cqq as a Lennard- Jones (LJ) sphere 



and the polythiophene chain as linked chain of Lennard- 
Jones beads with a bending potential between adjacent 
sites. The initial conditions were generated from a classi- 
cal molecular dynamics (MD) simulation of a single Cqo 
sphere embedded in cell with 25 polymer chains each with 
20 sites at constant T = 300K and volume as imposed 
by periodic boundary conditions. Upon photo-excitation, 
we assume that charge transfer occurs in the framework 
of a fixed geometry of the system producing a single Cq 
anion and a single positively charged chain located in 
close proximity to the the Cq anion. Since the chains 
are treated as a discrete set of beads, we use a bicubic 
spline to compute an interpolation function for 7(5) along 
the chain which is used for computing the curvature and 
Frenet-Serret frame along the polymer curve. We then 
assign local site charges by solving the Schrodinger equa- 
tion for a particle on a curve (Eq. [T]), adding to this a site 
energy determined by the Coulomb interaction between 
the sites on the chain and the anion, 

- ^L r + mn - i{s)i) - = Eip (3) 

where V(r) is the Coulomb interaction between the Cq 




(c)t = 200/ s (d)t = 300/5 

FIG. 2: Charge tunneling between kinks on a polymer chain. 
At t = the chain takes on a (+) charge and the Cqo takes 
on a (-1) charge. The color and magnitude of the sphere on 
the chain are indicative of the charge at a given site. At 
t = 20/s, the quantum state is clearly delocalized between 
kinks. At t = 30 fs the quantum state is localized at a new 
kink. 
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FIG. 3: Distance between the Cq and average + charge lo- 
cation vs. time. Inset figures show the polymer conformation 
and charge location at selected times. 

at location R and an arbitrary location along the chain 
7(s). We solved Eq. 4 numerically using a n-point 
discrete variable representation based upon Tchebychev 
polynomials. [TO] [IT] Site charges are given by \ip(si)\ 2 
evaluated at the center of each bead. We assume that 
the adiabatic approximation remains valid so that the 
interaction between the quantum and classical degrees of 
freedom can be described within the Hellmann-Feynman 
theorem. The information from the quantum calcula- 
tion was then passed to a molecular dynamics (MD) code 
(Tinker) [12] and for a single 1 fs MD step. This cycle was 
repeated for each dynamical step. Numerically, all this 
was accomplished using Mathematica (v. 8.0) [13] and the 
Tinker (v. 4. 2) molecular dynamics code. 

In Fig. [2ja-d) we show a sequence of frames from a 
representative simulation. At t = 0, shown in Fig. [2^i, 
we initiated the quantum dynamics by placing a non- 
mobile negative charge on the Cqq sphere and a mobile 
positive charge in the lowest energy eigenstate of a nearby 
polymer chain. The initial state is largely localized at a 
kink close to the Cqq. However, over the next 300 fs 
the charge tunnels from one kink to a nearby kink and 
becomes localized at the new kink. Similar transitions 
occur repeatedly over the course of a given simulation. 

In Fig. [3] we show the distance between the positive 
charge on the chain and the Cq over the course of the 
simulation. Inset figures give an indication of the geome- 
try of the chain and the distribution of charges along the 
chain at selected intervals. In particular, the dramatic 
jump at 15ps corresponds to the tunneling of the state 
along the chain towards sites closer to the Cq . 



In Fig. [4] we compare the effect of the various potential 
energy contributions in Eq. [3] at the same points of time 
highlighted in Fig. [3j The ground-state energy and wave 
function of each are also indicated in these figures. What 
is striking is that the curvature contribution (in yellow) 
is generally much less than the Coulomb energy between 
sites along the polymer chain and the Cqq . However, 
the presence of kinks on the chain give rise to deep traps 
along the chain that in many cases are stronger than the 
Coulomb interaction between the charge and the Cqq an- 
ion. In fact, in many cases, these kinks move the quan- 
tum state away from Cqq anion for long periods of time. 
This implies that in more realistic systems, dynamical 
fluctuations of the chains themselves may play a key role 
in the charge separation dynamics of bulk-heterojunction 
systems. 

Summary: The geometry of curves gives rise to an 
effective one-dimensional Schrdinger equation for a parti- 
cle bound to the curve in which the particle is attracted 
to regions of higher curvature. This purely a property 
of the geometry of curves and universal to all quasi-one 
dimensional quantum systems. We showed that when 
combined with a bending modulus, a particle bound to a 
linear curve will induce a spontaneous symmetry break- 
ing to bend the curve to optimize the curvature. 

We then applied this to develop a coarse-grained model 
of the charge-localization dynamics in a model of a 
Ceo^thiophene based photovoltaic material. Combining 
molecular dynamics of the polymers with quantum eval- 
uation of a mobile charge bound to a single polymer we 
show that kinks within the chain lead to localization and 
tunneling of the charge along the chain. We conclude 
that kinks and fluctuations of the chain may play an im- 
portant role in the charge fission and separation dynamics 
in polymer-based photovoltaic materials. 

These model simulations underscore the extent to 
which the dynamics and morphology of the polymer chain 
influences charge separation in polymer-fullerene based 
cells. While the current model does not account for quan- 
tum chemical effects, such as conjugation breaks due to 
orbital mis- alignment, chemical defects on the chain, or 
charge transfer between neighboring chains, such effects 
can be incorporated into the model in a systematic way 
and will serve the basis for subsequent investigations. 
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FIG. 4: Potential energy contributions and ground-state wave functions at various points along the simulation shown in Fig. 
3. Yellow = curvature potential, Red =Coulomb potential, Green = Total potential, Blue = wave function (scaled and off-set 
to ground-state energy). 



Physics) 62, 6296 (2000), URL |http : //link . aps . org/\ 



abstract/PRB/v62/p6296 

[2] J. J. Ladik, Fundamental World of Quantum Chemistry 
2, 271 (2003). 

[3] R. C. T. da Costa, Phys. Rev. A 23, 1982 (1981) 
URL jhttp : //link . aps . org/doi/10 . 1 103 /PhysRevA . 



23.1982 



[4] R. C. T. da Costa, Phys. Rev. A 25, 2893 (1982), 
URL jhttp : //link . aps . org/doi/10 . 1 103/Phys RevA ."| 
25.2893 

[5] J. J. M. Halls, K. Pichler, R. H. Friend, S. C. Moratti, and 
A. B Holmes, Applied Physics Letters 68, 3120 (19 96), 
URL |http : //link . aip . org/link/?APL/68/3120/l| 

[6] K. Harigaya, Physical Review B (Condensed Matter) 45, 
13676 (1992), U RL |http : //link . aps . org/abstract7| 
| PRB/v45/pl3676| 

[7] H. Imahori, Y. Sekiguchi, Y. Kashiwagi, T. Sato, 



Y. Araki, O. Ito, H. Yamada, and S. Fukuzumi, Chem- 
istry 10, 3184 (2004), ISSN 0947-6539 (Print). 
[8] N. S. Sariciftci, D. Braun, C. Zhang, V. I. Srdanov, A. J. 
Heeger, G. Stucky, and F. Wudl, Applied Physics Letters 
62, 585 (1 993), URL |http : //link . aip . org/link/?APL/| 
62/585711 

[9] J. Xue, S. Uchida, B. P. Rand, and S. R. Forrest, A pplied 
Physics Letters 84, 3013 (20 04), URL |http: //link, aip.] 
| org/link/?APL/84/3013/l| 

[10] J. C. Light, I. P. Hamilton, and J. V. Lill , The Journal 
of Chemical Phy sics 82, 1400 (1985), URL |http://link.| 
| aip . org/link/? JCP/82/14007T1 

[11] J. C. Light, in Time- Dependent Quantum Molecular Dy- 
namics (Plenum-Press, 1992). 

[12] J. W. Ponder, TINKER: Softw are Tools for Molecular 
Des ign, 4.2 ed (2004), URL |http : //dasher . wustl . edu/| 

| tinker/| 



5 

[13] Mathematica (Wolfram Research, Inc., Champaign, IL, 
2010), version 8. 



